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AN “AVERAGE” CONCEPT IN CALCULATING THE 
ELASTIC COEFFICIENTS FOR COMPOSITE MATERIALS 


Dumitru BOLCU’, Sabin RIZESCU’, Marcela URSACHE’, Cora MIREA* 


Rezumat. Aceasta lucrare propune..o..metoda noua, originala constand in 
considerarea unui concept de_,,medie’’ in calcularea.componentelor matricei 
rigiditatii in cazul barei compozite, cu doud faze, cand proprietatile elastice ale 
constituentilor sunt cunoscute. 


Abstract. This work proposes a new and original method consisting in considering 
an “average” concept in calculating the components of the rigidity matrix in case of 
two phases right composite bars, when the elastic properties of their constituents 
are known. 
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1. Introduction 


The composite materials are very useful in many industrial domains like aircraft 
and automotive. The elastic coefficients are determined mostly using experimental 
ways. A real challenge is the analytical calculus of these coefficients. 


Preliminaries 


Let’s consider a composite material made of some many distinct constituents 

(phases). Those phases could present different forms of anisotropy. 

Many authors [1], [2] take into account that under the action of certain external 

charges the material will accumulate for each phase: 

This way, the Sj, surface between the "1" phase and "k" phase (fig. 1) the 

continuity conditions are: 

- for the strain-stress status: 
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- for the deformation status: 
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deformations corresponding to three ortogonal directions n; ft ; 7 the way that: 


a 
n the normal unit vector to Six, from the “k” phase to the “1” phase. 


> > 
t and ¢ are unit vectors normal each other and they are situated in the 
> > 5 
tangent plane to Sin (tf x T =n ). 
The composite material can be assumed as being a continuum material. 
(i) 


Fig. 1 
An “average” concept for the calculus of the elastic coefficients 
for a two-phases composite material 

The basic concept of this workpaper [3], [4] is that if a composite has constituents 
having not too much different characteristics then a certain external charge will 
produce an unique deformation status as well as an unique strain-stress status that 
are parts of an unique constitutive equation. Both deformation status and strain- 
stress status have components with values that are to be calculated like some sort 


of averages of the values provided by fo} and je} that characterizes each 
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and every “i phase. This way, the material an be considered as being 
“homogenous” and it remains to determine its type of anisotropy. 
We’ll define: 


- the vector of average deformations {a} ; 
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k=1 V,. 


- the vector of average internal efforts { o} : 
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Where V is the volume of the composite, dV is the elementar volume 


dV =dx,dx,dx,, V, being the volume of the "k" phase [v = Vv, ' 
k=1 
Considering all these, the constitutive equation of the composite as a whole can be 


written as follows: 
(o} =[é] {2}: (5) 


or: {2} =[c}{o}; (6) 
where [E] and [C] are, respectively the rigidity matrix and the thickness matrix of 
the composite material. We’ll show that both matrices are symmetric, having their 
main diagonal different from zero, so they are inversable. 

In order to be more more specific, let’s consider a composite material 
presenting two homogenous and isotropic phases. In its non-deformed status the 
material is considered as having its own reference system Ox, x,x;. We note Sj2 


the separation surface between these two surfaces. We are to consider also a 
curvilinear coordinates system chosen the way that the in each and every point of 
Siz the unit vectors corresponding to this system are actually the same with the 
normal unit vector an other two rectangular unit vectors placed in the tangent 


> > > > > 
plane to Sj. We'll note these n; ¢ and 7c. If i, is the unit vector of Ox,; i, is 


3 
the unit vector of Ox,, and i, is unit vector of Ox,, so we have: 
> > > > > > > > > > > > 
n=,i,+ 8, i, + OR t =a,1,+ P,t, yONA: a0), 1, +2,1,4+06,i1,. J f7) 
Obviously, from (18) it results that the matrix of changing the unit vector basis 
from Ox,x,x, to the reference system attached to curvilinear coordinates will be: 


sl-le, B 3 |=|6, 2 B|: (8) 


Let’s consider: 


Onn On, Cnr O7; On O}3 
[o, |= Ont On O;, > [o]= O72 O72 O73 > (9) 
0, ‘s Oi, O07, O73 O73 033 


the (symmetric) matrices attached to stress-strain tensors in the local reference 
system defined by those curvilinear coordinates and in the Ox,x,x, reference 


system. 
Let’s also consider: 
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the (symmetric) matrices attached to.deformation tensors in those two reference 
system. Taking into account the tensorial character of the matrices from (9) and 
(10), the following relations does exist: 


lo}=[s]' lo, JIS): fe. I-[s] lo] Es1': (11) 
lel=[5} le, Is}: fe, F[sllellsT: (12) 


and, also: [s iv = [s li ? (13) 


[S] being an eigen matrix. Given all that we can write: 
to} =[a}to,} : to} =le] to} : (14) 
te} =[Blte,} 3 te, =[BI teh ; (15) 


and: 


where: 
\o} 19,3 F227 833,012, Guz, O23 } ? {o, } ACs O49 O29 Om One On j ? (16) 
\e} is 8095 8437 Vi ee {En } =F an Ty FE Ens Ens Eg ye (17) 


Concerning {o} and {e}, notations (16) and (17) are marking the order difference 
between the components 4 and 6 versus the usual notations and well known in 
literature. Concerning {o, } and {e,} we notice the same difference with respect 


to notation suggested by [o, | given by (9) and le, | given (10) which is to be: 


{o, }= {Fn , On , O,, , On , On , On, y > 1En es {Emm , Ent , on , en , Enr , Em y > dd 8) 

In order to find out the matrix of elastic coefficients we assume that, on the action 
of a certain external charge an internal status of strain-stress and deformation will 
appear in both phases. The column vectors attached to each status can be arranged 
in the curvilinear coordinates reference system as it follows: 


inl 
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where: 
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te, pyre Le, ey 272); 62)\, (22) 
and: 
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tT | as 
and, also: 
co, =f: 0 P60}; fo, f= fo? 02:02) (24) 


Considering all that we’ll come back to the two phases-composite and we 
can write for each phase: 


(eae ij 


leet eile ij 


where, materials being homogenous, we have noted the following matrices: 


and: 


i ee 3) UNO Le ae 
E; E; E; G; 
Dv. I Dv. I 7 
A,]=| NES 0. |: (ANG AS | Ie 0-431 =1,2 (27 
[4] mA ; [A ]=|0 : [c,] = i=1,2 (27) 
o Ned C7 49 0.<oAy 
G, | G, 


This way the constitutive equations for both phases take the following form: - for 
the first constituent: 


eV [A He, HIB Net}: 16, =1B: I Ye, 1c, Ho: }: (28) 


- for the second constituent: 


eV =TA. Ho, HIB tot}; te, J=1B. 1 fo, }+ IC Hof. (29) 


From (28) and (29) we conclude: 
\o ” \_ [A, li ie 2 [A, I’ [B, ea \ 
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and: 
io, ee : YA B,]\o; 7 
t,, J=1B re ‘ fe He} A) Le Me:}. G1) 
ae (30) an i, i is (29) and (20) we’ll have: 
ple lol fe; 
S ui , he = 
a [7;] lV; {e"} ' 
ele THES 7 
and: 
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be tie » 
where we used the following notations: 
RIB TAL > tla bedl-el| Hla }llstel| Slab) (36) 
I 0 0 0 0 0 
ih 1 i [o|= 0 i (37) 
00 1 0 0 0 
[D, |= [B, | [A, ies lz, Ic, [B, | [4,]" [B, ]; (38) 
[7,]=[4, Ie [uN=+[4, i [B, J: (39) 
[7, |=[4, ry lu, |--[A, Me [B, |: (40) 


Using (14) and (16) we can determine the column vectors for the strain-stress 


status and for the deformation status in the Ox, x, x, reference system: 


male aon eee lt ait a 
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Where, taking into account the concrete form for [a] and []: 
a; a; AQ; 
[M|=| BBS BB, |: 
20,8, 2a;8, a,B,+a,f, 
aa, aa; a 
[N]=| BiB: B,B; B; |: (43) 
a,P,+@,B, a,B;+a;B, 2a,f, 
[2a,6, 2a;6; a,6;+a,0, 
[P|= 2P,0, 2830; B,0;+ B30, |; 
EY * es 5,6, 
| HOMERS, @,0,+a,0, 2,0, 
[R]= £,0,+B,0, B)O;+B;0, 2,0, |. (44) 
| 6,6) 6,0; OF 
and: 
Ce oF 20a; 
El=| & BS | 2B,B. |; 
af, A,B; ,8;+a;f, 
2a,a, 2a,0; OF 
[A= 2f PB, 2f BP; 5 > (45) 
a)P,+a,B, a,B;+a;B, af, 
Ome C0, 4) 0, Og O40, 
[G]- B,0, P30; B20; + P30, |; 
Oo; We 20,0; 
G,0,+@,0, G,0; + Qh OF 
[H]= $0, + B,0, B)O;+P;0, Bo, |. (46) 
20,0, Lon, o; 


The column vector of medium deformations in the Ox, x, x, results from (3) and 
the column vector of the medium strain-stress status results from (4). Particularly, 
if {e, | and {o, } have constant components we obtain: 


Ti) Id i bly 
pa e] tr) fed te] 
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(47) 
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where: 
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(49) 


The constitutive equation for the composite as a whole in the Ox, x, x, reference 
system will take the following form: 


-\ [Al [B])s-), 
hier tel ts 
Putting in (50) \é \ given by (47) and \o } given by (48) we’ ll obtain: 


YM. bv, HoH} ¥ fe.) Ee J} st2l| Se. 1 tn) (51) 


k=1 
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The rigidity matrix contains both matrices [A] and [B] given by (51) and (52) and 
contains also both matrices [B’] and [C] given by (53) and (54). 
We obtain: 


a} [Sloe Weblo} Ste} tea) lS.) 


k=1 
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> > 
The unit vector n_ to the separation surface will have the Ox, direction, t will 


4+ 
have the Ox, direction and t will have the Ox, direction, so: 


n= i, > (a, =0; B, =0; 6, =1); ts i, > (a, =1; 8, =0;6, =0); t= : , 


(a, =0; B,=1;6,= ). 

100 0 0.0 00.0 1 0 0] 

a 1 Ao 0 ih Pp 0 a 1 0|; 

O° 4 00.0 Gaul 0 lo 01| 

1 0.0 00 0 00 0 1 0 0] 

or 1 br 0 sbi 0 2 1 0}: 

001 000 Onno jo o # 
With all that, we’ ll have 

Dll, HTM JP MEU Dale, |: in lo] 
= [7,.]=v, [47] TV; [A,}" = g| dL, | [z, =l0] 
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Seb ede», [2} 14, P' +» [8,140 =v: 
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Sik el-», fe, Hef) Bill fe. 1.1 [4.1 fe. = tl 


This way (51); (52); (53); (54) takes the form: 


[7] =[4llel: [ol Lally) [a lv} =[2' lel: le} =-L2'lly Ic]. 


Alle blab Lae 
slot’ bleb fal SPE EETT By, [A,)' BI: 
[a )-Wv) [ol =b, (2, [4,}' +y [B ie V+ [4b]: 
acres lic B14) Bs 
wileoHet a) Blkb Le eh | es tla} 
jv, [A,} +v, [A, ef [4,}'[B,]}+v, [4,7 [B, 


So: 
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Where v; and v2 are the volumic ratios for both phases. The matrices [Aj] and 
[Az] are symmetric. It’s easy to verify that: 

[al =[4]; [6 -[2]: [cf =[c]. (68) 
So, the rigidity matrix of the composite material is symmetric and that confirms 
once more the validity of the hypothesis adopted. 


So, we have: 
ome me E; (i-v3 )+v3 E; (L—v3)+2v,v, E, E,(I-v, v2). (69) 
' ' v, E, (elms a E, da DSS) 
. x VOB NE 5” )+v, v> BSW DS 
vf) =v) = KOM T 2. Oa? : i). (70) 
vy, E, (P= bv3/fy, F, j=" 
Gi}; =v, G,+v, G,; (71) 
bf Ke) le: Vy (1-v,)+v, V2 (7-v,)}- lv, E, (1+, )+v, E, (1+v,)] : (72) 
i . ip lt (1-v2)+v, E, =} 
G,G 
Gis = Cpa ae (73) 
1G, VIG, Wy 
— 2[v,v, 1—v,)+0, ¥, U-0,)P WE, (1+ oy) +v, E, (1+0,)] 


ne (7—v,)(1-v,) 3 25, (7-03 )+v3 ES (1—-v; )+ 2v,v, E, E, (I-o, v,)| 
Fo! (1+0))Z=2v,) | VD (7+, )(1-2v,) ; 


: (74) 
(I-v,)E, ( —D,)E} 
PE E 
(c) (c) 
V2) 1 023 
€ 7) 2 EO 0 0 va - E, Oi 
é Oo 
22 0 0 1 0 0 0 22 
Kal & Git O72|. (75) 
Ys1 0 0 0 at 0 go |}°% 
32 Gi; j O32 
Sis 0 0 0 0 oe 0 O33 
(c) (c) a ] 
Ss ee re) 0 6) —_— 
L £, E, ES” 


We have to specify that the "(c)" index refers to the entire composite. 


Finaly, we have, writing (75) in its usual form: 
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pio? pio? 
= nots 0 
E, Ei E, 
(c) (c) 
é a - ae 0 0 0 oO 
td Be ) EY ) E, 11 
E52 pio? vio? 1 O72 
13 23 0 0 0 . 
eae | E, E, EO O33 | (76) 
Y 32 I O32 
0 0 731 0 
Y31 G,; 03) 
Y21 0 0 0 0 ! O71 
Che 
vi 
0 0 0 0 
L Ge 
Conclusions 


The conclusion is that in this case the composite material can be considered as 
being homogenous and orthotropic having the planes Ox,x,; Ox,x,; Ox,x, as planes 


of elastic symmetry. 
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